Synchronization of mutually coupled chaotic lasers in the presence of a shutter 
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Two mutually coupled chaotic diode lasers exhibit stable isochronal synchronization in the pres- 
ence of self feedback. When the mutual communication between the lasers is discontinued by a 
shutter and the two uncoupled lasers are subject to self-feedback only, the desynchronization time 
is found to scale as A^t where > 1 and r corresponds to the optical distance between the 
lasers. Prior to synchronization, when the two lasers are uncorrelated and the shutter between them 
is opened, the synchronization time is found to be much shorter, though still proportional to r. 
As a consequence of these results, the synchronization is not significantly altered if the shutter is 
opend/closed faster than the desynchronization time. Experiments in which the coupling between 
two chaotic-synchronized diode lasers is modulated with an electro-optic shutter are found to be 
consistent with the results of numerical simulations. 
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Chaotic systems are characterized by an irregular mo- 
tion which is sensitive to initial conditions and tiny per- 
turbations. Nevertheless, two chaotic systems can syn- 
chronize their irregular motion when they are coupled 
0. When the coupling is switched off, any tiny pertur- 
bation drives the two trajectories apart. The separation 
is exponentially fast, and it is described by the largest 
Lyapunov exponents of a single system. 

In this Letter we show that the trajectory dynamics 
of coupled chaotic systems which also poses time-delayed 
self-feedback, is different. In a system with self-feedback, 
which has also been investigated in the context of secure 
communication with chaotic lasers [2(, the time scale for 
the separation of the trajectories is found to be much 
longer than the coupling time. On the other hand, when 
the coupling is switched on, rcsynchronization occurs on 
a faster time scale. We investigate this phenomenon 
numerically and show first experiments which support 
our numerical simulations. The demonstrated difference 
between de- and re-synchronization can be used to im- 
prove the security of public-channel communication with 
chaotic lasers 0- 

Semiconductor (diode) lasers subjected to delayed op- 
tical feedback are known to displays chaotic oscillations. 
Two coupled semiconductor lasers exhibit chaos synchro- 
nization. Different coupling setups such as unidirectional 
or mutual coupling and variations of the strength of the 
self and coupling feedback result in different synchroniza- 
tion states: the two lasers can synchronize in a leader- 
laggard or anticipated mode |3|, |^ , as well as in two dif- 
ferent synchronization states; achronal synchronization 
in which the lasers assume a fluctuating leading role, or 
isochronal synchronization where there is no time delay 
between the two lasers' chaotic signals 0, 0, El 13 • 

In this Letter we focus on a symmetric setup, the time 
delay between the lasers is denoted by r c and the time 
delay of the self- feedback is denoted by Td- In the event 



of t c = Td = t and for a wide range of the mutual cou- 
pling strength, a, and the strength of the self- feedback, 
k, the stationary solution is isochronal synchronization 
[3, 0- EJ • The quantity with which we measure the degree 
of synchronization between the two lasers is the time- 
dependent cross correlation, p, defined as 
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where I a and Ib are the time dependent intensities of 
lasers A and B and the summation is over times indicated 
by i. Isochronal synchronization is defined by the cross 
correlation, p, having a dominant peak at At = 0. 

We control the mutual coupling between the lasers by 
a shutter, located at a distance cr/2 from each one of the 
lasers, where c is the speed of light. When the shutter is 
open the two lasers are mutually coupled with strength a 
and with self-feedback k. When the shutter is closed the 
self- feedback, K e , is increased to a value of a + n so that 
the total feedback in the open/closed states remains a 
constant. This is required so as to prevent a sudden drop 
in the overall feedback, which would typically destroy the 
synchronization immediately |lfj|. 

The two quantities of interest in this letter are the 
desynchronization time, td, and the resynchronization 
time, t r . The desynchronization time is defined as the av- 
erage required time for the correlation to decay to Cdp(0) 
where Cd < 1. The time is measured from the moment 
the shutter is closed and p(0) is the average correlation in 
the isochronal phase. The resynchronization or recovery 
time is measured after the shutter has remained closed for 
a long period and the two chaotic lasers are uncorrelated, 
and K e = k + a. The shutter is then opened and the self- 
coupling strength is reduced to k. The resynchronization 
time is defined as the average time required, from the 
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shutter opening, for p to increase from zero to C r p(0), 
where C r is a constant < 1. 

To numerically simulate the system we use the Lang- 
Kobayashi equations that are known to describe a 
chaotic diode laser. The dynamics of laser A are given 
by coupled differential equations for the optical field, E, 
the time dependent optical phase, <E>, and the excited 
state population, n; 
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and likewise for laser B. The values and meaning of the 
parameters are those used in Ref . 0, IE IH] ■ The tunable 
parameters, both in the simulations and in the experi- 
ment, are k, a and p the pump parameter which is the 
ratio of the actual laser injection current to the threshold 
current. 

Figure 1 displays the calculated desynchronization 
time as a function of t for K = a = 50ns _1 , p(0) > 0.99, 
p=1.2, and Cd = 0.8 and 0.9. Each data point is av- 
eraged over 50 samples, and p(0) is measured by aver- 
aging a sliding window (sliding length is 1 ns) over a 
length r, while the solid lines are obtained by a linear 
fit, td — AdT + constant. The calculation shows that td 
scales linearly with r with a slope which increase as Cd 
decreases, and is near 7 and 9 for C d = 0.9 and 0.8, re- 
spectively. A similar linear scaling of the desynchroniza- 
tion time was obtained for all values of p in the range from 
1 to 1.5, where for a given Cd, t d decreases with p. The 
simulation assumes an ideal shutter with instantaneous 
closing and opening times, and also assumes a discontin- 
uous decrease of a to zero while k increases to k + a. We 
find that the linear scaling as well as the slope is robust 
to the following two experimentally necessitated pertur- 
bations: (a) a non-ideal shutter which closes gradually 
over a period of 10 — 20 nanoseconds; (b) an imperfectly 
closed shutter, allowing residual mutual coupling of a few 
percent of a while in the closed state. 

In the inset of Figure 1, the average p as a function 
of time for r = 50 ns is presented for a case where the 
shutter was closed at t=0. It is clear that the decay of 
p does not consist of a typical exponential decay. The 
striking result is that the correlation coefficient is almost 



a constant for a long initial period (first ~ 250 ns for 
the parameters of figure 1) and then crosses over to an 
exponential decay for very long times. Because the event 
of closing or opening the shutter takes a time r/2 to 
propagate to the lasers it is reasonable to expect that td I 
r/2. It is surprising, however, that td scales linearly with 
r with a prefactor which is significantly greater than 1. It 
is also not obvious from the simulation (in which p(0) > 
0.99), if such behavior can be observed in an experiment 
where p(0) < 0.9. 

The nearly constant value of p after the coupling be- 
tween the lasers is terminated calls for a theoretical ex- 
planation. Let us discuss the synchronization for the case 
k ~ a where both lasers are driven by an almost identi- 
cal delayed signal which is the sum of the self-feedback 
and the coupling beam. When the mutual coupling is 
switched off and replaced with stronger self-feedback, the 
system still feels its synchronized state for a period of 
length r, since the system is coupled to its history de- 
layed by a time r. The only difference caused by the 
closing of the shutter is that the lasers no longer commu- 
nicate with each other and each laser is coupled only to 
its own state. With time, a small difference in the driving 
signals develops. This small difference is amplified, since 
the system is chaotic, and this occurs stepwise in time in- 
tervals of length r. For each interval there is a constant 
distance between the two trajectories, which increases 
for the following interval. Only the envelope of these 
steps is described by the largest Lyapunov exponents, but 
not the dynamics itself. Hence desynchronization is very 
slow, and its time scale should be proportional to r. On 
the other hand, when the exchanged beam is switched 
on again, both lasers receive immediately an identical 
feedback signal (for k — a) and synchronize very fast, 
independent of r. We thus expect that the desynchro- 
nization time is insensitive to the values of n and a and 
should scale with r. In contrast the resynchronization 
time should be very sensitive to the difference, k — a. 

The experimental setup, which confirms many of these 
numerical perdictions, is shown schematically in Figure 
2. We use two semiconductor lasers emitting at 660 ran 
and operated close to their threshold currents. The tem- 
perature of each laser is stabilized to better than 0.01-ftT. 
The lasers are subjected to similar optical feedback and 
are mutually coupled by injecting a fraction of each one's 
output power to the other. A fast electro-optic modu- 
lator, with measured closing/opening time of 15 ns, is 
introduced in the middle of the coupling optical path to 
enable closing and reopening of the mutual coupling. The 
optical setup is designed so as to compensate the sudden 
drop in the overall feedback power when the shutter is 
closed and the mutual coupling feedback drops to near 
zero. We thus use the shutter as a polarization beam 
splitter which divides the output power of the laser into 
two parts: one used for the self-feedback and the other for 
the mutual coupling channels. The opening and closing 
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FIG. 1: td as a function of r for p=1.2, n — a — 50ns" 1 and 
C d = 0.9 (■) and C d = 0.8 (•). The lines are a result of 
a linear fit and the error bars obtained from 50 independent 
calculations for each point. Inset: The average p as a function 
of time for the same parameters as in the main figure and 
r = 50ns -1 . 




FIG. 2: Schematic of the two mutually coupled lasers. EOM - 
Electro-optic modulator; BS - Beam Splitter; PBS - Polarizing 
Beam Splitter; OD - Optical Isolator; NDF - Neutral Density 
Filter; PD - Photodetector. 

of the shutter merely changes the ratio of powers in the 
channels, but maintains the overall feedback power at a 
constant level. Without this precaution the sudden drop 
in feedback power would destroy synchronization imme- 
diately This setup, however, does not prevent the change 
in phase of the laser field when the shutter changes its 
state. This residual effect decreases the level of synchro- 
nization by a small amount (as can be seen in Figure 3). 
The shutter also does not close hermetically and the leak- 
age power to the mutual coupling channel in the closed 
state is measured to be ~ 7% of the shutter open value. 

The feedback strength of each laser is adjusted using 
a A/2 wave plate and an optical diode (see Figure (2J and 
is set to a value which leads to a reduction of about 5% 
in the solitary laser's threshold current 0. The lengths 
of the self-feedback and coupling optical paths are set to 
be equal to obtain stable isochronal synchronization 
Two sets of measurements are reported here, correspond- 
ing to two self-feedback optical paths with r = 12.7 ns 
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FIG. 3: Experimental results for p as function of t for p = 
1.08, 1.19 and r = 12.7, 23.6 ns . Each data point was 
averaged over 500 measurements. Inset: Data collapse, 1.7 -td 
(td) measured in ns as a function of p for r = 12.7 ns (23.6 
ns). 

and r = 23.6 ns. 

Two fast photodetectors (response time < 500 ps) are 
used to monitor the laser intensities which are simulta- 
neously recorded with a digital oscilloscope (500 MHz ,1 
GS/sec). The correlation coefficient, p, is calculated by 
dividing the intensity traces into 10 ns segments (each 
segment containing 10 points) and p is calculated be- 
tween matching segments and then averaged. 

The measured correlation coefficient, p, is shown in 
Figure 3, in which the shutter is closed at t=0. The 
coupling power decays to its closed level in about 15 ns, 
limited by the speed of the shutter. The observed decay 
time is only slightly shorter than the decay time obtained 
in simulations for the same r, and as in the simulations, 
p initially maintains a high and nearly constant value 
for 50 — 100 ns, which is much longer than r. The four 
data curves shown in the figure correspond to different 
experimental parameters, indicated by the value of p and 
to the two different values of r. 

Figure 3 shows that td increases with r and the inset 
of Figure 3, depicting td as a function p, indicates that td 
is a decreasing function of p. The inset of Figure 3 also 
demonstrates that the different decay curves all collapse 
to a single decay curve, independent of p, when scaled by 
a factor of 1.7 which is very close to the ratio of ri/ra = 
23.6/12.7^ 1.86. The numerical simulations for larger r 
also exhibit such data collapse when scaled by ~ T1/V2 
resulting in scaled decay curves which are independent of 
p|l3l|. This result and the linear scaling of td for a given 
p indicate 

t d (T, p) oc rg(p) (1) 

where g(p) is a function characteristic of the specific 
diode laser used. For small r finite size effects are ex- 
pected as a result of the positive constant in the linear 
scaling shown in Figure 1. Indeed, for r = 12.7 and 23.6 
ns, the numerical results indicate that the average ratio 
td(r\,p)/td(T2,p) ~ 1.75, which is in surprisingly good 
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agreement with the experimental result of 1.7. 

The simulations also indicate that a transition from 
the low frequency fluctuation (LFF) reg ime to the fully 
developed coherent collapse regime Il2j occurs at p ~ 
1.35, which is close to the experimental value of p (inset of 
Figure 3) where td becomes almost independent of nfl^. 
Though the decay time obtained from the simulations is 
longer that the decay time observed in the experiment, 
this is not surprising, since in the simulations the systems 
are initially correlated to a very high level (p(0) > 0.99) 
while in the experiments the initial correlation is p(0) < 
0.9 

We now turn to examine the scaling of the resynchro- 
nization time or the recovery time as a function of r. In 
the simulations we start with two uncoupled systems with 
self-feedback K e = n + a. When the shutter is opened at 
t = 0, K e is reduced discontinuously to k and the mutual 
coupling is set to a. For all examined cases, with k ^ a, 
our calculations indicate that the resynchronization time 
also scales linearly with r. This scaling is exemplified in 
Figure 4 for C r = 0.9, k = 60ns -1 and a = 40ns -1 . 

The inset of Figure 4 displays the resynchronization 
time as a function of k — a for a given t. It appears 
that this difference, rather than the coupling strength it- 
self, is what controls the resynchronization time, and t r 
scales almost linearly and symmetrically with k — a. For 
k = a, t r is very fast and simulations indicate that it is 
independent of r (limited by the fixed size of the sliding 
window) as expected. In all examined cases, the pref- 
actor of the linear scaling of the resynchronization time 
was found to be -C Ad, indicating t r < td- We also ob- 
served similar behavior, i. e. t r < td, in the experiment, 
though quantitative determination of the experimental 
resynchronization time is complicated by ringing in the 
high voltage electronics used to turn on the modulator 
and by the fact that the modulator response time is as 



long as 15 ns, which is comparable to t r . 

Although experimentally we cannot accurately mea- 
sure the resynchronization time, we show, in Figure 5, a 
demonstration of the persistence of the synchronization 
between the lasers upon repeated closing-opening opera- 
tions of the shutter. Shown is the typical behavior of p 
while the shutter is closed for ~ 40 ns and then reopened. 
The other parameters of the experiment are r = 12.6 ns 
and p ~ 1.08. The cross correlation coefficient p is not 
affected, by the closing/reopening of the shutter and the 
changing of n e and a, since td > 40 ns. 

The results reported above for re/de-synchronization 
times, which were also obtained recently for chaotic maps 
fl3| , demonstrate the possibility of establishing a reliable 
chaos based communication channel even while the com- 
munication between the lasers is interrupted by relatively 
long intervals. We expect that these effects will play an 
important role in advanced secure communications using 
mutually chaotic lasers. 
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